Artificiality of multifractal phase transitions 
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A multifractal phase transition is associated to a nonanalyticity in the generalised dimensions. 
We show that its occurrence is an artifact of the asymptotic scaling behaviour of integral moments 
and that it is not observed in an analysis based on differential n-point correlation densities. 
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Multifractal measures appear in a number of nonlinear physical phenomena like turbulence [|l||2j] , chaotic dynamical 
systems || and high-energetic multiparticle dynamics |Q, to name but a few. Due to the close analogy between 
the multifractal formalism and statistical thermodynamics any nonanalyticity in the generalised dimensions 

is interpreted as a multifractal phase transition pf]|-p^[. This behaviour has been discussed in the context of the 
above mentioned phenomena JToj , |l3| — |l5[| and is also denoted as the occurrence of strong intermittency. Note that 
conventionally generalised dimensions are extracted from the asymptotic scaling behaviour of moments, where the 
latter represent integrals over the fundamental correlation densities. We will now demonstrate that a multifractal 
phase transition is an artefact of the integral moment analysis and is not observed in an analysis based on differential 
correlation densities, where the true generalised dimensions of any order, characterising the underlying multiscale 
process, are revealed. 

As the textbook example of a multifractal process we consider a one-dimensional discrete random multiplicative 
cascade. It is associated with a binary tree structure, obtained by hierarchically partitioning the original interval of 
length lo = 1 into subintervals of size lj — 2~ J . The density e« 1 \ associated with a (j + regeneration interval 
characterised by the binary index k = (fei • • • fcj+i) with each k taking on possible values or 1, is multiplicatively 
linked to the density at the larger scale by 

0+i) =a [j+1) e U) (1) 

k\- ■ -kjkj + i "ki-'-kjkj + i k\---kj ' V / 



Independently for each branching, the left and right weights, qh = Qk^l-o an< ^ 1 R ~ Qk^-k-v are drawn from 
a probabilistic splitting function p(qL,qii) with support < qL,lR < oo. Without loss of generality we assume 
p(ql, qii) = p(qn, <7l) and set (q L ) = 1 as well as e (0) = 1. 

Nature does not allow an infinitely long multifractal scaling range; in fully developed turbulence, for example, it 
is restricted to r\ <C I <C L, where n and L represent the Kolmogorov and the integral length scales, respectively. 
Consequently, we restrict the random multiplicative cascade to a finite number J of cascade steps. The complete 
statistical information of the ensemble of generated cascade fields is then contained in the multivariate characteristic 
function 

Z[\^] = LJ t X^kAiu)) , (2) 

\ \fei,...,fej=0 J I 

from which the n-point correlation densities are derived 



(3) 

A(J)=0 



by taking appropriate derivatives with respect to the conjugate field variables X^. The multivariate characteristic 



function and the resulting n-point correlation densities have been calculated analytically in Refs. [[16 17|; see also Refs. 
For the extraction of generalised dimensions exponents so-called (box-) moments are considered, defined as 



M„ 



where the backward density 



has been resummed over the smallest scales from J to j. With (||) the (box-) moment (Q) can be understood as a 
(box-) integration over the n-point correlation density (g|): 

) = 2^71 E " E v^-kM^i •■•* ( , 1) '" e *^-***& 5 i-* ( r ) ) ' (6) 



k {1 2, ,...,k { ,'=0 k { "\ ,...,fc ( "'=0 
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yielding the explicit expressions: 



Mi(J,j) = l , 



(q L q R ) , 2 - (g 2 > - ( 9Lto ) 



q 2 ) 2-(«a) V 2 



(7) 



and, for arbitrary order, 



{p} £ri i= n ^ ' 

where {p} stands for all possible partitions of n,i = n with rij G {1, 2, . . . , n}; the coefficients are simple scale- 
independent functional of the splitting- function moments (?™ 1 <Zr 2 ) with < mi + TO2 < - While we will present a 
full technical understanding of the structure reflected in the expression M) at a later stage of this presentation, some 
intuitive understanding can already be derived: the backward density (|)|) can be rewritten as 

JJ,i) = Ji) . . . n U) 1 V F {J ~ j) 

kj+i kj=0 



Ji) 



where (1 + represents the resummed density of the subcascade with length J — j following the branching point 
(fci • • • kj). This resummed density need not be strictly equal to 1 and in general will fluctuate around 1. Hence, the 
backward density (||) need not be identical to the forward density e^..^ .- In view of (|TJ), the first factor (q")- 7 of the 
expression (||) then originates from ((sjjf^, ..&.)") while the remainder of the expression (§ is equal to ((1 + A^ J -^) n ). 

For the very special case of a conservative cascade, where with p(qL, qn) — p(Ql)^(ql + Qr. — 2) the sum of the left 
and right weight is conserved at every branching, the two coefficients a?\ of the second-order moment ([?]) become 
0{i,i} = (?i<?fl>/(2- (<? 2 )) = 1 and a {2 , } = (2- (<? 2 ) - (<7Lfe))/(2 - (<? 2 )) = 0, since (fete) = (q(2-q)). Similarly, all 
but one coefficients aj™, of the expression (||) vanish and the moment of order n becomes exactly M n (J,j) — (q™) 3 . 
Also in view of (^|) this outcome is intuitively clear since due to the conservative nature of the splitting function the 
resummed density 1 + A^ J ~^ is strictly equal to one. The moment M n (J,j) does not depend on the length J of the 
cascade and exhibits perfect scaling. The multifractal scaling exponents 

r(») = (10) 

are deduced by setting M n (J,j) — (/o/^') T ancl are related to the generalised dimensions D n by r(n) = (n — 1)(1 — 
AO- 

A factorised splitting function p(qL,<lR) — p(ql)p(qr) is a representative of non-conservative cascades, where the 
sum of the left and right weight is only globally conserved ((ql + <Z_r) = 2), but not locally [qi + qn ^ 2). As a 
consequence, the mixed splitting-function moments {q^qji) 7^ (9™ 1 (2 — g)™ 2 ) do not show the anticorrelation typical 
of conservative cascades. For this case coefficients Q^J\ are generally nonzero and the moments (||) do not show 

rigorous scaling. The resummed density 1 + &S J ~i> of the subcascade with length J — j now fluctuates around one 
and causes the deviations from rigorous scaling behaviour. 

For nonconservative cascades two subclasses of splitting functions have to be distinguished: in the case of so-called 
weak intermittency, the support of p(qL,lR) is restricted to < qL,QR < 2, so that the respective moments are 
restricted to (q n ) < 2"~ 1 , where n > 1 and where the extra 1/2 on the right hand side of this inequality comes from 
the requirement (q) = 1. This implies that, in the limit of an infinitely long cascade, J — > oo, only the true scaling 
term M n (J,j) ~ (q n )^ survives in the expression (||), so that asymptotically for j <C J the same multifractal scaling 
exponents ( [To| ) are extracted as from the corresponding conservative cascade. For the resummed density 1 + A( J ~^ 
of the subcascade, this implies that in this asymptotic scaling range its probability distribution has converged to a 
scale- independent fixed point; see also Refs. pp| , pT| . 

The second subclass of nonconservative splitting functions exhibits the phenomenon that has become known as 
multifractal phase transition or strong intermittency: once the support of p{qL,qR) allows values q^ and/or qn to 
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exceed 2, a critical order n cr it exists, so that (q n )/2 n ~ 1 < 1 for n < n cr it and (q n }/2 n ~ 1 > 1 for n > n a it- Then, 
again in the limit of a very long, but finite cascade and given that n > n cr ;t, the term corresponding to the partition 
{n, 0, . . . , 0} dominates the moment (H) of order n, which hence for j <C J scales as 



Mn>n„JJ,i<<J)«ag 0i ... i0} ^y 2 ^) ~ (I)"- 1 . (11) 



For the multifractal scaling exponents this implies 



rin< g «)/ln2 (n<n crlt ) 
v y \ n - 1 (n > ricrit) , v ' 

so that there is a discontinuity in the first derivative of r(n) with respect to the moment-order n at n cr it- For n < n cr j t 
the multifractal scaling exponents r(n) = ln(exp(n lng))/ In 2 may be interpreted as a free-energy-like function with 
moment order n as inverse temperature. - According to (p"l|), note that in the limit J — » oo of a non-physical, infinitely 
long cascade, moments with order larger than n cr jt would diverge. In view of (Q) this implies that the probability 
distribution of the subcascade-resummed density 1 + A^" 7-3 ' comes with an algebraic tail. 

For a nonconservative random multiplicative cascade with a factorised splitting function p(qL,qn) = p(ql)p(qr), 
where, for example, 




2-Koq 

is of log-normal type, the multifractal scaling exponents ( |l2] ) are found to be r(n < n cr j t ) = a 2 n(n — l)/(21n2) 
below the critical order n cr i t — 21n2/cr 2 , which defines the multifractal phase transition at r(n cr it) = n C rit — 1- In 
fully developed turbulence a good qualitative description of observed multiplier distributions in the surrogate energy 
dissipation field jH) and an acceptable intermittency exponent r(2) = 0.25 has been found for the parameter choice 
a = 0.42; for this value we find n cr it = 7.86. Note, that for the log-normal weight distribution the so-called Novikov 
rules p2]| do not apply since weights may exceed a value of 2 with nonzero probability. - The scale-dependence of 
the exact second-order moment expression (0) is illustrated in Fig. 1. As expected for j <g; J, M 2 {J,j) scales as (g 2 )- 7 
for a < <r cr it and as 2- J for a > <7 cr it, where a cl a = Vhi 2 corresponds to n cr it = 2. As j — * J noticeable deviations 
from this scaling behaviour occur. Directly at a = fT C rit, where (q 2 ) = 2, these finite size effects become so large that 
it becomes difficult to extract a scaling exponent asymptotically. Analogous findings hold for higher order moments 
M n (J,j). 

In the presence of a multifractal phase transition the information on the true multifractal scaling exponents ( fioj ) 
with n > n crit appears to get lost. But this is not the case! Note in this respect that, according to (§)-@, the 
moments M n ( J,j) are box-integrals over the n-point correlation density (^), the latter thus being more fundamental 
than the former. For demonstration, we pick the two-point correlation density 

pW{d) = = (l 2 ) J ~ d «to«R> + (!- (lLq R ))5 d , ) , (14) 

which is a function of the ultrametric distance d = J — j between the two bins K\ — {k\ ■ ■ ■ kjkj + \ ■ ■ ■ kj) and 
^2 = (hi ■ ■ ■ kjkj +1 • • ■ kj), where fej+i ^ k'j+i- K is depicted in Fig. 2 for a random multiplicative cascade with a 
factorised splitting function of log-normal type and reveals perfect scaling with the true multifractal scaling exponents. 
For the two-point correlation density and, in general, for all n-point correlation densities no multifractal phase 
transition occurs. Then, why does a multifractal phase transition occur once based on integral moments and not 
on correlation densities? 

The answer to this question will be found in an inconspicuous property of the correlation densities. For demon- 
stration we discuss this again only for second order. First, we realize that the second-order moment (^) can be cast 
in the form 

M 2 (J,j) = ^- (^pW(d = 0) + E 2 d ~ V 2l (rf)) ■ (15) 
Next, by adding and again subtracting something for d = 0, the expression (O) is rewritten as 
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P m (d) = ( q 2 y 



VI l <in) + ( 2 _ ^ - (<7;,<7/e, 



(<7 2 > - (qLQR/ , i..; 



2-<<7 2 > 



<<zT&o 



= p u (d) + 2 _ (^ <9 ) <5do • (16) 

For d ^ the modified two-point correlation density fP^ (d) is identical to the original two-point correlation density 
pi 2 \d). The difference comes for d — 0, where 

pP](d = 0) = (q 2 ) 

pPl(d=l) 2-( g 2 ) 1 j 

contrary to pl 2 l(<i = 0)//9l 2 l(e? = 1) = (q 2 ) / (qLqii) ■ If we were only to substitute the modified two-point correlation 
density into (|l5|), then this difference insures perfect moment scaling with the true multifractal scaling exponents as 
we arrive at the first term of the expression (R). Substitution of only the second term appearing in (16), which is 



proportional to Sdo, produces the second term in ([?]), which is proportional to {{q 2 ) /2) J 2 J and reflects a trivial scaling 
with exponent r(2)| tr iviai = 2 — 1. The appearance of the second term in (^) and (|l6|), respectively, is a consequence 
of the missing anticorrelation between weights gz, and qn in the splitting function; it vanishes only for a conservative 
splitting function p(qL,qii) = p(qL)S(qL + qs. - 2) because then {qhqR} = 2— (q 2 ). 

It is worthwhile to elaborate in more detail on the occurrence of a multifractal phase transition from the viewpoint of 
Eq. (|l6|). For demonstration we pick again a factorised splitting function of log-normal type (|l3|), where (q 2 ) = exp(cr 2 ) 
and (<7Lte) = 1- As a increases monotonically from to a CI a = Vhi2, the ratio (g 2 )/(2 — (q 2 )) appearing in 
( |l7| ) increases from 1 to +oo; then, as a is further increased, it changes sign and increases from — oo at a = cr^ it 
monotonically to —1 as a — > oo. Consequently, as the (d = 0) and (d = 1) elements of the modified two-point 
correlation density contribute as the sum p' 2 '(d = 0) + p^{d = 1) to the second-order moment ([l5]), they enhance 
each other for a < cr C rit, but more or less cancel each other for a > <r cr jt- Hence, for the former case the modified 
two-point correlation density dominates the moment scaling whereas for the latter case the (S-function-like correction 



in (16) becomes dominant. 

We conclude: when generalised dimensions are determined via integral moments, one is likely to encounter artificial 
multifractal phase transitions, which are not a property of the underlying strongly intermittent nonconservative 
cascade process, but artifacts of small-scale resummation. The fundamental n-point correlation densities and, should 
the underlying process only be resolvable at an intermediate scale, so-called density correlators 



<(#>) 


""■-(«"») 


> 


{<&>) 


"">■■■((*")'""> 



/>l"«l>-">™»nj _ \ \ / \ / / /ION 

^ "((^' ) ) mi )-((^T") ' 

avoid such contributions and are therefore a better choice in estimating generalised dimensions. 
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FIG. 1. Scale-dependence of the second-order moment M^J-, j) resulting from a finite discrete random multiplicative cascade 
with a factorised log-normal splitting function. The dashed straight lines represent the asymptotic scaling M2(J,j) ~ (q 2 ) j 
and M2(J,j) ~ 2 j for a < cr cr i t = %/rn 2 = 0.833 and a > a cr i t , respectively. 
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